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ABSTRACT 



We report on H-band, ground-based observations of a transit of the hot Neptune GJ 436b. Once combined to achieve sampling 
equivalent to archived observations taken with Spitzer, our measurements reach comparable precision levels. We analyze both sets of 
observations in a consistent way, and measure the rate of orbital inclination change to be of 0.02±0.04°in the time span between the 
two observations (253.8 d, corresponding to 0.03±0.05°yr~' if extrapolated). This rate allows us to put limits on the relative inclination 
between the two planets by performing simulations of planetary systems, including a second planet, GJ 436c, whose presence has 
been recently suggested (Ribas et al. 2008). The allowed inclinations for a 5M e super-Earth GJ 436c in a 5.2 d orbit are within ~7°of 
the orbit of GJ 436b; for larger differences the observed inclination change can be reproduced only during short sections (<50%) of 
the orbital evolution of the system. The measured times of three transit centers of the system do not show any departure from linear 
ephemeris, a result that is only reproduced in < 1 % of the simulated orbits. Put together, these results argue against the proposed planet 
candidate GJ 436c. 
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1. Introduction 

Among the known transiting exoplanets, GJ 436b is by far the 
smallest and least massive exoplanet, being to date the only tran- 
siting hot Neptune. With a mass of 23 M®, it orbits a M star with 
0.44 M Q , in a slightly eccentric orbit (e=0.15), with a period of 
2.64 d. The exoplanet GJ 436b was discovered by a radial veloc- 
ity survey (Bu tler et al.L l2004). and in the same work the authors 
repo rted a non-detection of photometric transits. Three years 
later, iGillon et all (|2007b ) announced the discovery of transits of 
this object, which provided the first (and currently the only) mea- 
surement of the radius and density of a Neptune-sized transiting 
exoplanet. Its radius is currently reported to be of around 4.22 R® 
(Torres, 2007), and according to models o f icy and rocky plan- 
ets dSeager et al.ll2007HFortney et al.ll2007h, it is thought that the 



plane t has an H-He envelope (Gil lon et aill2007bl iDeming et al.l 
2007). Several authors have signaled that the significant non- 
zero eccentricity is an effect of a further companion in the sys- 
tem that excites the eccentricity, otherwise t he system should be 
circularized in rela tively short timescales (Ma ness et al.ll2007t 
IDeming et ap2007l) . 

The low i nclination of t he or bit (86.3°), the non-detection 
of transits by iButler et all (|2004|) . and an analysis of a low- 
significance peak in the residuals of the radial velocity data led 
Riba s et alj(l2008l) to suggest the presence of a 5 M s super-Earth 
(GJ 436c) orbiting with a period of 5.2 d. An important argument 
of this hypothesis is the change in the inclination of GJ 436b in- 
duced by the gravitational interactions with GJ 436c. Under the 
assumption that there was a ~0.3 deg variation in the orbital in- 
clination between 2004 and 2007 -which would have made the 
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* Based on observations taken with the Telescopio Carlos Sanchez 
(TCS) of the Observatorio del Teide, operated by the Instituto de 
Astroffsica de Canarias. 



transit non-detection bv lButler et aL ( 20041) and the d etection by 
IGillon et all d2007bl) compatible- iRibas et ail d2008l) suggested 
that the transits observed in 2008 should last ~2 minutes longer, 
due to a rate of inclination change of roughly 0. 1 °yr _1 . 

In this work, we report ground-based H-band observations 
of a transit of GJ 436b observed in March 2008, and perform 
a rean alysis of the Spitzer 8 //m data reduced by IGillon et all 
(2007a). We use different techniques aimed at measuring the rate 
of inclination change on GJ 436b, and provide constraints on the 
inclination and masses of the proposed planet GJ 436c. 



2. Observations and data reduction 

Data were collected at the 1.52 m Telescopio Carlos Sanchez 
(TCS) telescope, operated by the Instituto de Astroffsica de 
Canarias at the Observatorio del Teide, during the night of 
March 8, 2008. We used a H filter and the detector CAIN-lfl 
a NICMOS 3 technology chip, which employs a 256x256 array 
of HgCdTe elements sensible to the range 1-2.5 microns. The se- 
lected W widefield optical configuration of the detector provides 
an image scale of 1 "/pixel. 

Due to the high brightness of the star, we severely defocussed 
the image in order not to saturate the star with a 1.5 s expo- 
sure time. This resulted in ring-shaped images of the star with 
an outer diameter of ~26" and an inner diameter of ~1" . The 
flux of the star was thus spread into ~500 pixels, which helps 
to minimize the effects of bad pixels, tracking, and atmospheric 
seeing changes. We selected a part of the detector with few hot 
or dead pixels and, to achieve a good stability of the system, 
we employed no dithering pattern. To account for hot pixels in- 
side the target's PSF, we moved the telescope every 1-2 hours 
to a close zone of the sky with no stars and took 50 exposures. 



http://www.iac.es/telescopes/cain/cain.html 
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Fig. 1. Normalized light curve of GJ 436 on the night of March 8, 
2008. The dispersion of the data out of transit phases is 0.0023, 
with a median time sampling of 2.35 s. 



Pixels with a low or zero sensitivity (dead pixels) were located 
in images of the dome under a slight artificial illumination. They 
represent ~3% of the detector. A Fowler readout mode was em- 
ployed, performing 8 readouts per image and rejecting 2 of them. 
Taking into account the observing overheads of typically 1 s, we 
obtained a total of ~8000 exposures in a single night. 

To extract the flux of the star, we first interpolated the dead 
pixels we found as explained above. We used the average images 
of each set of 50 background exposures to subtract the hot pixel 
contribution. The flux was summed inside a circular aperture of 
20" in radius and centered on the star, and the background level 
was estimated from an annulus of 22" and 34" with the same 
center. The integrated flux of the star in a single 1.5 s expo- 
sure was ~2.6xl0 7 e~. We restricted our analysis to the phases 
close to the primary transit of GJ 436b. The star exhibited a long 
term variation of only ~ 1 % during the night, which was well 
corrected by a parabolic fit to the parts outside the transit. In 
two moments of the night (around HJD-2454534 = 0.5736 and 
0.603) we observed a fast decrease and increase of the flux, in 
a 20-30 s timescale, probably caused by slight vignetting of the 
dome. We discarded these points for the rest of the analysis. The 
final normalized light curve that was used for the analysis de- 
scribed in the next section, containing 3600 points, is plotted in 
Fig-Hi The remarkably low noise level attained can be attributed 
to both the smooth variations of the atmospheric transparency 
in the H filter, as well as to the decision not to use a dither- 
ing pattern, which could have emphasized the detector inhomo- 
geneities. 



3. Analysis 

In order to evaluate the rate of orbital inclination change for 
GJ 436b, we modeled the TCS and the 8/mi Spitzer data in a 
consistent way. Both observations of transits are separated by 
253.8 d. For th e transit modeling, we used the formalism of 
iGimenez] J2006I) . which is also valid for eccentric orbits, such 
as the case of GJ 436b. The best fit was found by a minimiza- 
tion of the y 2 u sing the commonly used AMOEBA algorithm 
dPress et al.ll 19921) . To estimate the error of each of the fitted pa- 
rameters, we performed 100 fits in different sets of light curves, 
which were constructed by subtracting the best-fit model, shift- 
ing circularly the vectors of the residuals and their errors by 
a random number, and again adding the best-fit model. With 
this bootstrapping procedure, we allow for the possible low- 
frequency structure of the residuals due to uncorrected system- 
atic effects. 
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Fig. 2. TCS H-band phased light curve of GJ 436 and best-fit 
model (green line), and the residuals from the best-fit model 
(bottom). The bin size corresponds to ~27 s, and the 1-sigma er- 
ror bars have been estimated from the dispersion of the points in- 
side each bin. The standard deviation of the residuals is 0.00077. 



3.1. TCS data 

We grouped the TCS data in bins with a width of ~0. 00012 in 
phase, corresponding to ~27 s, which is approximately the same 
time sampling as the Spitzer data set. The errors were estimated 
as the standard dispersion of the points inside each bin, divided 
by the square root of the number of points inside each bin. The 
solution that provides the best fit is presented in Table [TJ The 
fitted limb darkening coefficients are in good agreement with the 
values reported by IClaretl d2000l) for a star with T c// =3500 K, 
logg=4.5. 

3.2. Spitzer data 

One transit of GJ 436b has been observe d with Spitzer at 8 jum 
dDeming et all 120071; iGillon et alj|2007al) . We analyzed the set 
of data reduced by IGillon et al] d2007af) . obtaining the results 
in Table [T] We note that we obtain a different inclination than 
that reported bv lTorres et al.l (2008), which is lower by 0.2°. By 
modeling the transit with a zero eccentricity, we obtain an in- 
clination in perfect agreement with the inclination reported by 
Torres et al. (2008), and ~0.4°higher than the solution proposed 
by Gillon et al. (2007a). In Fig. [3] we plot the three different fits, 
which shows the importance of properly including the eccentric- 
ity of the orbit in the transit modeling. Due to the low-frequency 
noise in the parts of total eclipse, we obtained unrealistic values 
of limb-darkening if we also fitted for thes e parameters. We thu s 
fixed these parameters to those reported bv lGillon et al.l (|2007a). 



3.3. Comparison of both data sets 

One of the predicted effe cts of a planet with similar character- 
istics as that proposed by Ribas et al. (2008) is the variation of 
the orbital inclination (and thus, the total duration of the tran- 
sits) of GJ 436b. Measuring the absolute orbital inclination us- 
ing transits is a challenging task, due to the effects of uncer- 
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Table 1. Parameters of the GJ 436 fits, and associated one-sigma 
errors. 
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Fig. 3. The region of ingress/egress of the Spitzer data of 
Gillo n et all (|2007a), and three different models: 1) Our best-fit 
model, with an inclination of 86.54°and taking into account the 
eccentric orbit (solid green line). 2) Our best-fit model, for a cir- 
cular orbit (dott ed blue line) . The i nclination is 86.3°, in perfect 
agreement with lTorres et al.l(l2008b . 3) A model with an inclina- 
tion of 86.3°and taking into account the eccentric orbit (dashed 
red line). Note that for eccentric orbits, the times from the transit 
center to the transit ingress and egress are slightly different, as a 
consequence of the change in orbital velocity during the transit, 
making the lines appear doubled in this plot. 
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tainties on the limb darkening coefficients and the possible ef- 
fects of active regions on the star. A simpler approach consists 
in measuring relative orbital inclination changes of the same ob- 
ject by the more precise -and less affected by the uncertainties 
on the limb darkening coefficients- determination of the total 
transit duration. As a first, simple method, we evaluated the to- 
tal duration of the two transits and their transit centers with a 
fit to a trapezoidal function, where the fitted parameters were 
the transit center, duration, depth, and ingress time. The fit was 
perfo rmed with a Levenberg-Marquardt algorithm dPress et alj 
1992). The duration of the transit observed with Spitzer is es- 
timated, with this approach, to be 58.7+0.7 m, while the TCS 
best-fitted trapezoid lasts 59.2±0.8 m. The centers of the tran- 
sits were used to determine the ephemeris of Table [TJ The er- 
ror bars in the timing of the transit centers were estimated from 
the dispersion of the bootstrap fits described in the previous sec- 
tion, and thus include the influence of uncorrected systematics 
in both sets of data. Alternatively, we can also compute the total 
transit duration from the value of the fitted 9\ ; in this case, it is 
necessary to perform the conversion between the circular orbital 
phase used intrinsically in the IGimenez (2006) formulation and 
the observed orbital phase. Using this method, we obtain a dura- 
tion of 60.5±1.0 m and 61.0+0.7 m for Spitzer and TCS transits, 
respectively. We attribute the slight difference between the ab- 
solute values using both methods to the effect of the limb dark- 
ening, considered to be zero when performing trapezoidal fits. 
The relative difference in total duration AD of the transits using 
both methods is in perfect agreement. This difference in dura- 
tion can be converted to a difference in orbital inclination Ai () t, s 
of 0.02±0.04°between the two observations, or 0.03+0.05°yr~' 
if we assume the Ai b S behaves linearly in these timescales. 



4. Discussion 

In order to put physical limitations on a second planet from the 
observational constraint of the previous section, we assembled 
numerical models of two-body planetary systems. We defined 
the initial osculatin g Kepl e rian o rbital element^ of the GJ 436b 
from the values of Torres] d2007l) and chose the orbital plane of 
this planet as the reference plane of the system. The elements 
of GJ 436 c were defined according to the values proposed by 
Ribas et al. (2008). The parameters for the planet c were chosen 
randomly in the intervals resumed in TableffJ the angular param- 
eters not present in this table were chosen randomly for the two 
planets. The evolut ion of the syste m was calculated using the 
RADAU integrator (Everhart ll985t) . which ensures a proper re- 
production of the close encounters between the planets. We per- 
formed a total of 5 000 simulations, following the evolution of 
the system for 1400 yr, which are ~2xl0 5 orbits of GJ 436b, and 
~10 5 orbit of the outer planet. The orbital parameters of GJ 436b 
were recorded every 253.81 days, which is the difference in time 
between TCS and Spitzer observations, to be able to perform 
comparisons between the Ai obtained from the simulations and 
the observed Ai t, s . For each of the simulations containing N time 
stamps, we calculated the fraction / = N{Ai < Ai b S )/N of time 
stamps in which the change of the inclination of GJ 436b was 
within 1 <x of the measured Ai t, s . We explored the variation of 
/ as a function of the mass and inclination of the second planet. 
We created a regular grid in this parameter space, and computed 
the mean value of / in each cell. The result is represented in 
Fig. |4] As expected from orbital interactions, massive planets 



2 The observed Keplerian orbital elements at a particular epoch. See, 
for instance. lRovUl988h 
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Table 2. Mass of the star GJ 436, initial osculating Keplerian 
elements for GJ 436b and ranges for GJ 436c. 



GJ436 








M, 


0.452 M Q 






GJ 436b 




GJ 436c 




in h 


23.17 M m 


m c 


1-6 M e 


(ib 


0.02872 AU 


a c 


0.044-0.046 AU 


e b 


0.15 


e c 


0-0.2 


k 


0° 


ic 


0-20° 


U) b 


343° 


w c 


222-308° 




2 3 4 

Mass of second planet (M E ) 

Fig. 4. The fraction of time that the simulated planetary systems 
verified the Ai f, s restrictions summarized in Table [TJ as a func- 
tion of the mass of the second planet and its orbital inclination. 
The horizontal lines mark the observed inclination of GJ 436b 
(solid line) and the lower limit to the inclination of GJ 436c to 
produce transits (dashed line). 



induce relatively large inclination changes unless their inclina- 
tion is close to the inclination of GJ 436b. Thus, the change in 
inclination sampled by only two points will be less than Ai b s 
only in a very small fraction /. A 5 M @ planet has a value of 
/ bigger than 0.5 only when its inclination is within ~7°of that 
of GJ 436b, and planets as less massive as 3 M ffi produce com- 
parable values of / if their inclination difference is around 15°. 



planet, this fraction was of 0.5%, rising to 5% for the 2M ffi case. 
Hence, there is a very low probability that the non detection of 
TTV in three epochs is in agreement with the existence of a 5M e 
companion in a period close to the 2: 1 resonance. 

However, the low-measured value of Ai b s , compatible with 
zero within 1-cr, and the absence of significant TTV both 
strongly argue against the proposed 4.8M ffi planet in a 5.2 d or- 
bit. Further measurements of transit epochs could serve to re- 
strict the few remaining systems that produce apparent TTVs of 
the order of the observed values. 

Additionally, we have shown that ground-based observations 
with a 1.5 m class telescope in the near infrared can attain preci- 
sions that are of the order of those obtained with Spitzer at higher 
wavelengths. These kind of observations are specially suited for 
transit modeling and TTV studies. 
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Furthermore, the suggested period for GJ 436c is close 
to a 2:1 mean-motion resonance, what would induce tran- 
sit time variations (TTV) of the order of several minutes 
dHolman & Murravll200H lAgol et al.l l2005). For two particular 
cases (5M ffi and 2M ffi ), and assuming a difference in orbital in- 
clinations of 5°, we calculated the times of minimum projected 
separation between the centers of the star and the planet GJ 436b, 
to track the evolution of the TTV for 10 yr. In order to compare 
with the observed values of center times (Table [TJ, we randomly 
chose a start epoch Eq, and calculated the standard deviation 
of the TTVs at Eq,E^2\, and E+^f, (there are 21 orbits between 
the discovery transit of Gillon et al. 2007b and the transit ob- 
served with Spitzer, and 96 between this last transit and the tran- 
sit observed with TCS). We repeated the process 10000 times, 
and calculated the fraction in which this value was smaller than 
the standard deviation of the observed TTV (12 s). For a 5M e 
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